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Abstract 

We provide new information concerning the pseudospectra of the 
complex harmonic oscihator. Our analysis illustrates two different 
techniques for getting resolvent norm estimates. The first uses the 
JWKB method and extends for this particular potential some results 
obtained recently by E.B. Davies. The second relies on the fact that 
the bounded holomorphic semigroup generated by the complex har- 
monic oscillator is of Hilbert-Schmidt type in a maximal angular re- 
gion. In order to show this last property, we deduce a non-self-adjoint 
version of the classical Mehler's formula. 
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1 Introduction 



The harmonic oscillator is known to be ubiquitous in many branches of quan- 
tum theory. Unlike most quantum mechanical problems, it can be solved 
explicitly in terms of elementary functions (see section 2), so it provides an 
excellent illustration for the general principles of the quantum-theoretical for- 
malism. It also plays a key role in quantum electrodynamics and quantum 
field theory, and has been a focus of great attention among the mathematical- 
physics community for many years. 

Although it might seem very difficult to say something new about the 
harmonic oscillator, in a couple of recent papers ([^ 0]) E. B. Davies pro- 
vides a framework which opens the possibility of obtaining new information 
concerning the spectral theory of this and other Schrodinger operators. In 
the present paper we extend some results given in and and introduce 
some new ideas that will be useful in studying both spectral and stability 
properties of this operator when it has a complex coupling constant. 

We define the complex harmonic oscillator to be the operator 

HJix) :=-^f(x) + cx'fix) 

acting on L^(]R) with Dirichlet boundary conditions, where c G C is such 
that Re(c) > and Im(c) > 0. Our main interest is to investigate the 
resolvent norm of Hf. inside its numerical range. To this end, if we let the 
e-pseudospectra of He be 

Spec,(ife) :=Spec(i/,) VJ{zeC: \\{H,- z)-^\\ > ^} 

for all e > 0, our interest is in obtaining information concerning the shape of 
these sets for small e. 

We adopt the notation and basic ideas about pseudospectra recently de- 
veloped by E. B. Davies, L. N. Trefethen and others. For general results 
and more examples in the present spirit, we refer to ||, |], |12| and to the 
bibliography there. 

The results we obtain here extend in two ways those given by E. B. Davies 
in where he shows that for any positive constant h independent of?] > 0, 

||(i/e-r/[6 + c])-i -.oo 

as ^ oo. On the one hand in section 3 (theorems ^) we demonstrate that 
for all 6 > and 1/3 < p < 3 fixed. 
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as 1] —>■ oo. On the other hand in section 5 we prove (as a consequence of 



theorem 12) that for all 6 > there exist a constant > such that 



lim \\{Hc-[r] + ib]y \\ < Mb 

and 

lim \\{H,-[cri-ib])-^\\ < M^. 

r/— >oo 

The study of pseudospectra provides important information about the 
stability of operator He- By virtue of the identity 

Spec,(ife) =\J{SpeciH, + A) : P|| < ^ 

(for a proof see |11|]), knowing the size of these sets for e close to zero, 
allow us to obtain very precise information about the stability under small 
perturbation of the eigenvalues of He, 0. It will turn out (sections 3 and 
5) that high energy eigenvalues are far more unstable than the first excited 
states (see also 1T|, [§])• 

The method we use in section 3 is analogous to the techniques in p, |^ 
and involves the construction of a continuous family of approximate eigen- 
states for He by means of JWKB analysis (see 0] for an introduction to this 
topic and some applications). We should mention that such a procedure is 
examined for more general potentials and from the numerical point of view 
ini. 

In section 4 we deduce a non-self-adjoint version of the classical Mehler's 
formula, i.e. we construct an explicit formula for the heat kernel of —He 
(theorem |^. This formula allows us to show that the bounded holomorphic 
semigroup generated by —H^ is compact in a maximal angular domain. The 
results of section 4 provides a new approach to get resolvent norm estimates 
for He- Based on these estimates, in section 5 we obtain a set that encloses 
SpeCg (He) for small enough e (corollary [TT] and theorem |T2p. This confirms 



the numerical evidence given in H] about the shape of such sets. 



2 Definitions and Notation 

We will suppose in this paper that the parameter c G C satisfies Re(c) > 
and unless explicitly stated we do not impose conditions under Im(c). 

We assume that He acts in L^(R) with Dirichlet boundary conditions as 
follows. Take the closed m-sectorial quadratic form 

Qc{f,9)-= / f'{x)g'{x)dx + c / x'^f{x)g{x)dx 
Jr Jr 
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for all f,g e W^^^{R) n {/ G L\R) : / x^\f{x)\'^dx < oo}. Then is 
defined to be the m-sectorial operator associated to Qc via the Friedrichs 
representation theorem (see [0). 

The subspaces C^(M) and S{M.) (the Schwartz space) are form cores for 
He- Since 

C^(M) C S{R) C W^'2(M) n {/ G V{R) : j x^\f{x)\^dx < oo}, 

it is enough to check the desired property for C^(M). By |TD|, theorem 1.21, 
p. 317] a subspace is a core for an m-sectorial form, if and only if it is a 
core for its real part. Notice that the real part of Qc is the non-negative 
quadratic form Re(Qc) = <5Re(c), thus using for example p, theorem 1.13] we 
can deduce that C^(M) is a core for Re((5c) and therefore that it is also a 
core for as needed. 

If c is positive, is self-adjoint, but if Im(c) 7^ 0, is not even a normal 
operator. Moreover (see 0) there does not exist an invertible operator U 
such that UHcU~^ is normal. 

If we put Xn '■= c^/^(2n + 1), the spectrum of He is 

Spec (He) ={A„ : n = 0, 1, ...}. 

It consists entirely of eigenvalues of multiplicity one and it is routine to check 
that if Hn is the n^'^-Hermite polynomial, then 

vi/„(x) := c'/'H^{c'/^x)e-^'"'"^"^/^ 

are the eigenf unctions of He, so that 

if,^„, = A„^„ (n = 0, 1, ...). 

For 0<a,/3<7r/2, we shall denote the angular sector 

S{-a,(3) := {z e C : -a < arg(z) < /?}. 

We will also put S{a) := S{—a,a). It is clear from the definition that 
the numerical range of the operator He is contained in the angular sector 
5(0,arg(c)). 

Proposition 1 For all Re(c) > 0, 

Num(/fc) = {h + eta G C : ti, > 0, tits > 1/4}. 
//Im(c) ^ 0, then 

Spec (He) C Int(Num(i7c)). 
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Proof Since «S(M) is a core for Qc, 

mm{H,) ={Qc{f) : 11/11 = 1, / e 5(M)}. 
By the Heisenberg inequality, we know that for any / e S{W) 

(^Jjf'{x)\'dx^ (^l^x'\f{x)\'dx^>l\\f\\\ 

Therefore clearly 

Num (H^) CW :^ + ct2 e C : ti,t2 > 0, t]t2 > 1/4} . 

Let us check the reverse inclusion. For this we need to find test functions 
/ e 5(M) such that 

ge(/)e W = {t + ^eC : i>o}. 

For t > 0, let 

ft{x) := e-*^' G 

Using elementary properties of the Gamma function we can calculate 

/ \f^{x)\^dx = 2t27rV2(2i)-3/2^ 

J^xV,{x)\'dx = ^(2^)-^/^ 
Combining this three equalities, we obtain for alH > 

yjn{x)\'dx^ ^l^x'\Mx)\'dx^ ^ 

Therefore 

and so every point in d{W) is in Num (He) . Since both sets are convex we 
have a fortiori 

W C Num (He) . 
For the second part, suppose that Im(c) 7^ 0. Since 

Spec (iJj C Num (H^) , 
5 



it is enough to show that 

d{W) n Spec (i/e) = 0. 

The eigenvalues of He he on the hne c^^^r, r > 0. The boundary d(W) only 
intersect this line at the point 



c 



1 1/2 



Clearly |-2o| < |c^^^|, therefore the desired property follows from the fact that 
zq is never an eigenvalue of He- I 

From the basic theory (see for instance 0), we know that Spec^ (He) 
contains the ^-neighbourhood of Spec (He) and it is contained in the e-neigh- 
bourhood of Num (He) . We will show later (corollary |lT]) that if Im(c) > 0, 
there exist an E > depending on c such that 

Spec, (He) C Num (He) 

for all e < E. Observe that if c is real, since He is self-adjoint, Spec^, {He) 
is actually equal to the e-neighbourhood of Spec {H^ , we stress that this 
property is false in general for non-self-adjoint operators. 

Proposition 2 For lm(c) > fixed, the resolvent norm of H^ is symmetric 
with respect to the axis c^^'^r, r G M. As a consequence for all e > the 
e-pseudospectrum of He is also symmetric with respect to this axis. 

Proof For all a > let 

TJ{x) := a'/^f{ax) {x G M), 

the operator is isometric in L^(R) and such that = T^-i. Putting 
arg(c) := "d, since 

T|^|i/4-f/'e''»^c|-V4 = |c| ^^"^Hc, 

it becomes evident that we can assume without lost of generality c = e*'^ for 
0<^< n/2. 

For simplicity we rewrite the operator 

He = P^ + cQ^ 
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where Pf{x) := if'{x) and Qf{x) = xf{x) (respectively the quantum me- 
chanical observables of momentum and position). Thus, applying Fourier 
transform we obtain 

\\{P' + cQ'-z)-'\\ = \\{Q^ + cP^-z)-'\\ 

= ||(p2 + e-*V-e-*''^)-i 

= ||[(p2 + e-*V-e-*''2)-Tll 
= \\{P' + cQ^-e'''z)-'\\ 

for all z ^ Spec {He) , which is precisely our claim. ■ 



3 High Energy Eigenvalues 

In this section we show that if the coupling constant c is such that Im(c) > 0, 
and Zr^ G Num (He) parameterized by > is such that 

\zri — {bq + cq^) \ ^0 (as ?7 ^ oo) 

for some 6 > and 1/3 < p < 3, then 

lim ||(ifc - Zr,Y^\\ = oo. 

Such a result will be a consequence of theorem ||. Notice that by proposition 
H it is enough to assume that 1/3 < p < 1. 

Our first aim is to obtain test functions frj G C^(R), parameterized by 
r] > 0, such that if 

1/2-7/2 I 2 7 I 2 
Zr, = tcT] ' " + a 7]' + ca 7], 

where a > and 1 < 7 < 3 are fixed constants independent of r/ > 0, 

We follow a similar procedure as the one given originally in [^. Let 

$(x) := e-*(^) {x G M), 

where the polynomial 

^(xo + s) := ^is + ^2sV2 + ^3^73 {s G R), 
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centered in xq '■— ar]^^^, has coefficients 



ZC 'v/9 / 

= ——T] ^' [1 + cq 



l-7\ 



Clearly 

H^^{xq + s) = + Zri)^{xQ + s) 

where 

p{s) = CiS + C2S^ + CgS^ + C4S^ (s G M) (2) 

has coefficients 



Cl = 


— V 




a 


C2 = 


0, 




^2 


C3 = 
















C4 = 





Let us establish some properties of the function $ above defined. A 
straightforward calculation implies that 

|$(xo + = exp [-/32(?7)s' - (^z{il)s^] , 

for 

/32(r7) := \m{c)7]^'^'^'/^ > 0, 
/33(r?) := 



Im(c)(l + 2Re(c)r7^-^) ^_^/2 ^ q_ 
So; 



Taking derivatives of \^{xq + s)|2 with respect to s, allows us to conclude 
that this function has a local maximum at s = , and a local minimum at 

2f32{v) ^ 26^77^/2 

3/?3(r?) ~ 1 + 2Re(c)r;i-^ ^ °° 

as 7/ — > CX3. 

The required can be defined truncating the function $ as follows. It is 
routine to define a C°°(M) compact support function 

9{x) = 





if 


\x — Xq I 


< r]^o 




if 


a; — 





for 60 := 7/6, such that there exist constants gi > and q2 > independent 
of X, 1], a or 7, with the property 

\g'ix)\ < q^r,-'\ 
\g"{x)\ < q2V-''\ 



for any a; G M. Then, 

/^(x) := g{x)(l>{x) (x G M). 

The next two lemmas point out some properties of which can be em- 
ployed to demonstrate (0). The constants for A; = 1, 2, . . . below are real 
and independent of 77, but can possibly depend on a or 7. 

Lemma 3 For as above and 1 < 7 < 3, there exist positive constants ai, 
02 and E^, independent of rj, such that 

for all rj > E^. 

Proof Since 60 < 1/2, 

2(33(ri)ri^H^ 2(1 + 2Re(c)r]^-^)ri^H ^ 
lim = lim -7T = 0, 

uniformly for < t < 2. Therefore, there exist E^ > such that for any 

< ff^M'ft^ < >h^, (3) 

for any < t < 2. 

For the first inequality: If 7/ > E'^, we have 

m? = I \g{xo + s)\^\^xo + s)\^ds 
Jm. 

Jo 
Jo 

> / exp[-2(32{vW^'t'^W"dt 
Jo 

> / exp[-2Im(c)M>-(^-^)/^rfM 

> ai7?-(i-^)/^ 
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Observe that we are using the fact that (5o + (1 — 7)/4 = 1/4 — 7/I2 > 0. 
The second inequahty is similar. ■ 

We now estimate the numerator in the left hand side of ([l|). By (0) we 
have 

\\HJ,,-zM < \\2g'ixo + s)^'{xo + s)\\ + \\g"{xo + s)^xo + s)\\ 

4 

+ ^||cfcs'=<l>(xo + s)^(xo + s)||. (4) 

k=l 

Lemma 4 Let > as in lemma^. There exist positive constants a^, 05 
and ag, independent of t] > 0, such that 

\\2g'{xo + s)<^'{xo + s)f < 0477^^/6 exp[-^^r/(=^-^)/6] 

\\g"(xo + sMxo + s)r < a577-^/^exp[-^^r/(3-^)/^] 
\\s''^xo + s)gixo + s)f < aev^^-'^^"'^'^^^ 
for all 7] > and k = 1, 3 or 4. 

Proof We use similar arguments as the one provided in the proof of lemma 
I- 

Let 

n:= {s eR : ri^° < \s\ < 2ri^°}. 
For all s G fi, we have 



2|2 



a 



Then for all rj > E^, by the conditions imposed on g{x) above and by 
we have 

\\2g'{xo + s)<^>'{xo + s)r = 4 [ \g'{xo + s)\'\^\xo + s)\'mxo + s)\'ds 
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Ji 2 
< a4r/^^/^xp[-^^r/(3-7)/6]. 

The second estimate is similar so that we can find without difficulty. 
Finally for r] > and A; = 1, 3 or 4, using again 



s^''mxo + s)\''ds 

Jo 

Jo 2 

Notice that (|T]) can be easily obtained from lemma 0, lemma ^ and equa- 
tion (I). Using such result, we can achieve the following theorem. 

Theorem 5 Let He be the complex harmonic oscillator such that Re(c) > 
and Im(c) > 0. If 

Zrj = brj + erf 

where b > and 1/3 < p < 3 are constants independent of rj > 0, then 

lim ||(ifc - = oo. 

Proof Assume l/3<p<l. If we put the unitary operator for a > as 
in the proof of proposition recall that for arg(c) := 'd 

T|g|i/4ii/gii5T|^|-i/4 = |c| ^^^Hc- 

For all r > and /5 > let 

Wr,,r := I3rj + (3re''^T]P. 

Since 

\\{H,.,-w,,i,\)-'\\ = \\{\c\-'/'H,-w,,\,\)-'\\ 

= |c|V2||(if,-|c|VX|c|)-^||, 
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using (1) putting c? = \c\^^'^P and thinking of |c| = r, it is clear that 

hm \\{HqH> — Wr,r)^'^\\ = c>o 
jj— »oo ' 

for all r > and P > 0. 

Now if we put P = |c|^/^6 and r = \c\/b, we obtain 

\\{H,-z,r'\\ = \\{\c\'/'H,.. - z,)-'\\ 

as r/ ^ OO. The proof can be now completed by proposition 0. ■ 

This theorem has some consequences for the £-pseudospectra of He- It 
implies that for e > fixed, any curve Zrf G Num (He) parameterized by 
T] > such that 

\zrj — {hq + crf) \ — (as — > oo) 

for some 6 > and 1/3 < p < 3, will eventually be inside Spec^ {He) and 
it will stay there as ^ oo. In particular this shows that high energy 
eigenvalues are increasingly unstable under small perturbations. 



4 Non-Self- Adjoint Mehler's Formula 

In this section we show that the bounded holomorphic semigroup of contrac- 
tions generated by —H^ is compact in a maximal angular sector. In order 
to deduce this property for the semigroup, we obtain an explicit formula 
for the heat kernel of He as in Q and show directly that this kernel is of 
Hilbert-Schmidt type. 

We say that for a, /3 > 0, a parameterized family of bounded operators Tj. 
in the Banach space is a hounded holomorphic semigroup of contractions 
in the sector 5'(— if and only if: 

1. T^i+r2 = Tt-^Tt-^ for all G S{—a,(3). 

2. ||T,|| < 1 for all r G 

3. Tt- is a holomorphic family of operators in r G S{—a,(3). 

4. For all / G i3 and £ > 0, 

limT./ = / 

for r inside S{—a + j3 — e). 
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It follows directly from the definition that, for —a < ^ < (3 fixed, Tei^j for 
t > is a Co one-parameter semigroup in the standard sense, [Q. The gener- 
ator of Tt- is, by definition, the infinitesimal generator of the one-parameter 
semigroup Tt where t > 0. 

For convenience, we shall put 



Theorem 6 Let He the complex harmonic oscillator for Iie{c) > as defined 
above, then: 

1. For c fixed, —He is the generator of a hounded holomorphic semigroup 
of contractions e~^'''^ , with parameter r in the open sector Int(5'c). 

2. If Im(c) > 0, then iHc and — e*('^/^~^''s(^))ffc are also generators of a 
one-parameter semigroup such that e"^"'^ is strong continuous for all 

res,. 

3. //Im(c) < 0, then —iHc and ~e~^^'^^'^~^^^^^^^'' He are also generators of a 
one-parameter semigroup such that e^^"'^ is strongly continuous for all 

res,. 

4- For r > fixed, e'^""^ is a holomorphic family of bounded operators 
parameterized by c for all Re(c) > 0. 

Proof Properties 1 is deduced without difficulty from [|1^, theorem 1.24, 
p.492] and @, theorem 2.24]. 

Property 2 and 3: if Im(c) > 0, operator iH, is maximal dissipative 
operators, therefore is the generator of a one-parameter semigroup (the same 
argument works for the other three cases). From analyticity, follows strong 
continuity for all r G Int(5'c) and strong continuity in the edges can be 
checked using [||, corollary 3.18]. 

Observe that if Im(c) = 0, Spec {±iHc) is purely imaginary, so we cannot 
apply the above argument. 

Property 4: by the way in which we define its domain, the operator H, 
is a holomorphic family of type (B) for Re(c) > 0. Using the fact that 
holomorphic families of this type are locally m-sectorial, this property can 
be demonstrated by analogy to fT^, theorem 2.6, p. 500]. ■ 
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arg(c), f) \ {0} if Im(c) < 0. 



-arg(c))\{0} if Im(c) > 
) if Im(c) = 
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For T E Sc let the coefficients: 

A := exp[-2c^/V; 
Wi := 



W2 



W3 



c^/^(l + A^) 
2(1 -A2) ' 



(1 - A2) 



For all X, y G M, put 

Kc{t, X, y) := wi exp [w^xy - W2{x^ + y^)] , 
and define the integral operator 

^c,r/(a;) := / Kc{T,x,y)f{y)dy, (5) 



for all / G L^(R) for which this formula makes sense. 

Fixing r > 0, Kc{T,x,y) is holomorphic for Re(c) > 0, and fixing c it is 
continuous for t E Sc and analytic in its interior. Observe that for c real, 
Kc{t, X, y) has poles if r runs along the complex axis; this situation is avoided 
by the definition we have made of the angular region Sc- 

As pointed out in modifications of the classical argument (see 
show that 

e-^^^ = Ac,r (6) 

for all c > and r > 0. Using analytic continuation this equality can be 
extended to complex c and r. 

Theorem 7 (Non-self-adjoint Mehler's Formula) Let Kc{t, x,y) be de- 
fined as above. Then for all Re(c) > and r G Sc, 

e-^-/(a;) = / Kdr, x, y)f{y)dy (/ G 
Jr 

The rest of this section is devoted to proving this theorem. First we 
establish some local bounds for \Kc{T,x,y)\, this allows us to show that Ac^r 
is a holomorphic family of Hilbert-Schmidt operators in both parameters and 
then we obtain the desired equality by analytic continuation. 

Lemma 8 For all tq G Sc, there exists a neighbourhood Tq G Vq C C and 
real constants a\, a2 and such that for all t eVq{~\ Sc: 

\Kc{t, x,y)\ < ai exp [a^xy — a2{x^ + y^)] (allx, y G M). 

The constants satisfy: ai, 02 > and 2a2 ± 03 > 0, locally uniformly in c 
and T. 
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Proof Clearly 



\Kc{T,x,y)\ = |wi| exp[Re(w3)xy - Re{w2){x^ + y^)]. 

Since Wi are continuous in r and c, it is enough to show that Re('u;2) > and 
Re(2w2 ± W3) > 0, for all Re(c) > and r G Sc- 

Let := arg(c^/^) and 'd^- := arg(r). Without loss of generality we will 
assume that Im(c) > and |c^/^| = 2. Define the Moebius transforms by 

M^{z) {ze C). 

Then 
and 

2«;2 ± W3 = c^/^ii^ = 2e^'^^M±(A). 
1 =F A 

Notice that M4. maps the interior of the unitary disk into the open right 
half plane, and multiplying by e*'^'^ rotates about the origin by i&c- Every disk 
centered in the origin, of radius ri < 1, changes under M+ into a disk D2 of 
radius 

2ri 
1- r( 

centered at 

C2 := 1 > 0. 

If ri is close to 0, r2 is small, C2 is close to 1 and the rotation of D2 about 
the origin by 'dc remains in the right half plane. If ri is close to 1, r2 is big, 
and rotating by i}c can send some points of D2 to the left half plane. In spite 
of this possibility, notice that points in D2 which are to the right of C2 do 
not cross to the left half plane when rotated. It is elementary to find the 
maximum radius ri which allows e^'^''D2 stay in the right half plane. We call, 
the disk with center the origin having this radius, the critical disk. 

Now if we consider -^c and i9r fixed, describes a spiral with radius 
decreasing exponentially as |r| increases, starting in 1 and coiling about the 
origin. We leave to the reader, to check that all these spirals cross the critical 
disk with real part large enough to guarantee that after the mapping by M+ 
and the rotation by 'dc, the resulting curve stays in the right half plane. This 
gives Re(w2) > and a similar analysis gives Re(2w2 ± W3) > 0. ■ 
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Lemma 9 For all Re(c) > and r G Sc, the operator Ac^r as defined in ^ 
is of Hilhert- Schmidt type and 

1. For c fixed, Ac^r is norm continuous for t E Sc and holomorphic for 
T e Int(^c). 

2. For T > fixed, A^^j. is holomorphic family of bounded operators for 
Re(c) > 0. 

Proof Lemma ^ implies that for all tq G Sc, there exist a neighbourhood 
To G Vq C C and a constant M < oo such that: 

J J \KciT,x,y)\'^dxdy < M, 

for all r G Vo n Sc, where the constant M depends locally uniformly on 
c. Therefore is Hilbert-Schmidt operator and its norm depends locally 
uniformly on c and r. 

By the dominated convergence theorem, it is elementary to check that 

\\Kc{T,;-)-Kc{To,;-)h^O 

as r — > Tq. This provides norm continuity of Ac^r in r G for c fixed. Finally 
the analyticity of the family of operator in both variables, can be deduced 
without difficulty by differentiating under the integral sign. ■ 

By theorem ^ and lemma we know that the families of bounded opera- 
tors e~^'='^ and Ac^r are both holomorphic in each parameter. They coincide 
when those parameters are real. If we fix r > 0, by analytic continuation in 
c, equation (||) is also true for all Re(c) > 0. Now fixing c with Re(c) > 0, 
by analytic continuation in r, the same equation is true for all r G Int (S'c). 
Finally, since both families are strongly continuous in r at the edges of Sc, 
we have proved theorem ^ as requested. 

Lemma ^ is one of the crucial points for the analysis of pseudospectra that 
we intend to carry out in the next section. In particular, we use strongly the 
fact that e~^'='^ is a compact operator for all r G Sc- 

5 Pseudospectra 

Our aim in this section is to employ the compactness of the bounded holomor- 
phic semigroup generated by —He (see previous section) to obtain estimates 
on the resolvent norm of He inside S'(0, arg(c)). The technique is based on 
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the spectral radius formula for the semigroup and a convenient Jordan de- 
composition of the operators involved. For simplicity we will assume from 
now on that Im(c) > 0, but with a few corrections, the same results and 
techniques apply to the case Im(c) < 0. 

We establish first a formula for the spectral radius of one-parameter semi- 
groups; a proof of this fact can be found in theorem 1.22]. If e~^* is a Co 
one-parameter semigroup, the limit 

a := lim log ||e~"^*|| 

always exists with — oo < a < oo, and the spectral radius 

rad (e-^*) := max {|A| : A G Spec (e"^*) } = e"* 
for all t > 0. This implies in particular that for all a < —a 

lim e"*||e-^*|| = 0. (7) 

t— >oo 

Applying this estimate to e'^"'^ for r on the edges of Sc, we can show 
that the resolvent norm of is uniformly bounded in lines parallel and close 
enough to the edges of 5(0, arg(c)). 

Theorem 10 For fixed Im(c) > 0, consider the complex parameters Ziow = 
rj + ie and z^^pp = cij] — ie)/\c\, with rj > and e > 0. Then, for < d < 
Im(Ao) there exists a constant M^d < oo, independent of rj and e, such that 

\\{H,-zxo^)-^\\ <Me,d 

and 

for all rj > and < e < d. 

Proof By theorem |^ and lemma the bounded holomorphic semigroup 
e~^'''^ is compact for all r in the maximal sector Sc- Therefore, 0], 

Spec (e"^^") = {0} U {e-^"" : n = 0, 1, . . .}. 

By proposition § we can just concentrate on the lower edge. Putting r = —it 
for t > 0, we obtain 

rad(e*^^*) = e'^'^^^o)*. 

Fix < d < Im(Ao); by formula (0) with T = —iHc, a = — Im(Ao) and 
a = {d — a)/2, there exists t^ > such that 
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\\ds 




POO 


e"rfs + 






J Ta 




POO 


e^'ds + 


/ e^'^-'^^ds 








POO 


e^'ds + 


/ g-(Im(Ao)- 






. m 





for all t > ta- Then 

\\{H,-{r] + te))-'\\ = \\{iv-e-tH,)-'\ 

< 

< 

< 

< 

< W,d 

This result provides information about the shape of the e-pseudospectra 
of operator He- The result below will be extended in theorem 0. 

Corollary 11 For all < 6 < 1 there exists an e > such that 

Spec, (H,) C S{0, arg(c)) + 6c^^^. 

As a consequence of this corollary together with proposition |I|, we obtain 

Spec, (H,) C Num (H,) 

for e small enough. 

The situation for the other eigenvalues, requires more care and the argu- 
ment involves a Jordan decomposition of the problem. Let Qn the spectral 
projector of operator He associated to the eigenvalue A„, i.e. 

Qnf:=7^ [ {z-H,)-'fdz 



2iTi 



In 



where 7„ is a smooth curve whose interior just contains eigenvalue Xn- Ob- 
serve that Qn are not orthogonal projections in L^(]R). 
For m = 0, 1, . . . put 

m 

Pm ■ ^ ^ Qn- 

n=0 

It is clear that Pm is the spectral projector with rank m associated to eigen- 
values Ao, . . . , Am- Notice that L^(M) can be decomposed as the direct sum 
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of the closed subspaces Ran {Qn) for n < m and Ran (/ — Pm) , in the sense 
that the subspaces are hnearly independent and 

L2(M) = Ran (Qo) + • • • + Ran (Q^) + Ran (/ - . 

It is easy to show that each of the m + 2 subspaces above is invariant 
under the semigroup e~^'''^ for all r G Sc- Even more, the generator of the 
bounded holomorphic semigroups of contractions 



e-^^1 



and 



Ran(Q,i) 



|Ran(/-P™) , 

are respectively -i/c|Ran(Q„) and -i/c|Ran(7-p„) • Notice as well that, by 
compactness of the restriction of the semigroup, we have 



Spec (e-'^nRMi-Pm)) = W U {e-'-^}^=, 



m+l- 



In order to extend corollary |TT| beyond the first eigenvalue, observe that 
for all z ^ Spec (He) 

m 

[H, - Z)-' = Y^iHc - z)-^Qn + [H, - Z)-\I - Pm). 
n=0 

Then if we call Km '■= 1 + X]™=o HQ" 11' ^ straightforward calculation allows 
us to obtain the following estimate: for all z ^ Spec (He) 

(m 
'^\\iHc\K^n(^Q„) - z)-^\\ + \\{Hc\K^n(I-P^) -z)-^\ 
n=0 

(9) 

With the help of this inequality, we can achieve the theorem below. 

Theorem 12 Let He the complex harmonic oscillator such that Re(c) > 
and Im(c) > 0. For all < 6 < 1 and m = 0,1, . . . there exists an e > such 
that 

m 

Spec, (H,) C [5(0,arg(c)) + (A„+i ~ S)] U [j{z e C : \z - Xn\ < S}. 

n=0 

Proof To estimate the first sum of norms in (^, since Ran (Q^) is the 
one dimensional subspace generated by the n*'* eigenvector \E'„, the operator 



19 



-f^c|Ran(Q„) acts Oil this subspace as the operator of multiplication by An. 
Then for all z ^ X^-, 

||(-f^c|Ran(Q„) " ^)~"^|| 



I An — ^1 

To estimate the last resolvent norm in (j^), we use equation (H), applying 
an analogous of theorem |10] to -f^c|Ran(/-p„) instead of H^. Notice that now 
the first eigenvalue is Am+i. ■ 



In the notation of Aslanyan and Davies, [| 

m 

1 + ^k(A, 



where /t(An) is the index of instability of the eigenvalue An- Based on the 
results in [§], as n increases, the indices of instability ^(An) grow faster than 
any power of n. This is reflected in the above theorem in the fact that as m 
gets bigger the e we must choose gets exponentially smaller. 

This theorem confirms the numerical calculations made by Davies in 0, 
where he uses a computer package to find level curves for the resolvent norm 
of a discretization of the operator H(.. 

Going beyond the scope of theorem our conjecture is as follows. Let 



< p < 1/3, m = 0, 1, . . . and < 5 < 1 be fixed. Let 6„i p > such that 
there exist E > {] (possibly depending on m or p) verifying 

hk,pE + cE^' = Xm. 

Put 

Zrj ■■= bm,pr] + erf [t] > 0) 

and 

^m,p ■= {\zn\e'^ eC : 7]> E, arg(z^) < 1^ < arg(cz;;/|c|)} . 



Conjecture 13 There exists an e > sueh that 

m 

Spec, {H,) C VLrr,,p U{j{z eC : \z-Xn\< 5}. 

n=0 

This would be a substantial improvement of the results provided in this 
paper. The case p = is precisely theorem Because of theorem ^ the 
statement is false for l/3<p<lsoin this sense the constraint < p < 1/3 
is optimal. 
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